Single-file diffusion and kinetics of template assisted assembly of colloids 
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We report computer simulation studies of the kinetics of ordering of a two dimensional system 
of particles on a template with a one dimensional periodic pattern. In equilibrium one obtains a 
re-entrant liquid-solid-liquid phase transition as the strength of the substrate potential is varied. 
We show that domains of crystalline order grow as ~ i 1//z , with z ~ 4 with a possible cross-over to 
z ~ 2 at late times. We argue that the i 1 ^ 4 law originates from single-file motion and annihilation 
of defect pairs of opposite topological charge along channels created by the template. 

PACS numbers: 74.25. Qt,61.43.Sj,83.80.Hj,05. 65. +b 



Template assisted ordering is an useful method of 
preparing large, micro-arrays of functionalized nano- 
particles for a variety of technological applications A'~ 
3 . While there has been considerable amount of work 
on the structural and functional aspects of such arrays, 
their formation kinetics is relatively unknown in spite of 
obvious theoretical and technological interest. An im- 
portant aspect of structure formation in such systems, 
as we show in this Letter, is the possibility of coupling 
between low dimensional transport and phase ordering 
in a higher dimension. Specifically, we report on the 
ordering kinetics of a two dimensional (2d) system of 
colloidal particles placed in a one dimensional (Id) pe- 
riodic potential arising either from a template pattern 
etched on a substrate[3] or produced using crossed laser 
beams [4--9J. The periodic potential produces narrow 
parallel channels along which the motion of particles is 
severely constrained, since particles cannot go past one 
another without climbing the crests of the applied poten- 
tial. Such constrained motion of particles or single file 
diffusion (SFD), which also occur within narrow pores 
and channels [10] . eg. in zeolites [IT], carbon nanotubes 
or ion channels in cells has garnered a fair amount of 
attention involving experiments [12], theory [13] and com- 
puter simulations [14] . One of the characteristics of SFD 
is the development of long-ranged correlations in particle 
trajectories such that the mean squared particle displace- 
ment at late times t behaves as t 1 / 2 rather than the linear 
law for normal (Fickian) diffusion 13J. 

The freezing of a 2d solid, with lattice parameter 
a, in a commensurate, Id periodic potential^ V{y) = 
Vq cos(27ry/d) of wavelength d = y3a/2 has been the fo- 
cus of rather intensive recent research. Initial mean-field 
theories of laser induced freezing[T5] predicted re-entrant 
behavior where one obtains, at a suitably high density, 
first a freezing transition to a 2d triangular solid, as the 
amplitude of the laser field is increased, followed by re- 
melting at still higher laser intensities into a strongly 
modulated liquid phase. While the first transition is 
driven by a co-operative effect of density modulations in- 
duced by the field, re-melting occurs due to a reduction 
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FIG. 1: (color-online) Phase diagram from NAT Langevin dy- 
namics simulations of 10 4 particles interacting via the pair re- 
pulsive WCA potential in an external periodic potential in the 
density p and amplitude (3Vq plane at temperature T = 0.5. 
Blue/dark-gray open circles mark the state points for a liquid 
(L), solid (S) and re-entrant liquid phases. The (red/light- 
gray) arrows show two quench protocols where the system 
was first equilibrated either in the L phase (Qi) or in the re- 
entrant liquid phase (Q2) and subsequently quenched to S. 
Upper inset: Triangular lattice of lattice parameter a with 
the position of the crests of the external potential of wave- 
length d — v3a/2 marked by horizontal parallel lines. Lower 
inset: order parameter <f) of the system (see text) as a function 
of f3Vo plotted using open triangles for p — 0.87. The blue 
(dark-gray) dotted line is drawn at the cut-off value of cj> i.e. 
Sbd> = 0.35. 



of dimensionality from 2 to 1 as particles becomes in- 
creasingly confined along the lines of maxima in the laser 
intensity. This phenomenon was subsequently verified 
in simulations [5] and experiments [5]. Finally, re-entrant 
laser induced freezing in two dimensions was explained 
using a defect mediated mechanism in [7]. The validity 
of the defect mediated approach has been verified in de- 
tail by direct computer simulations [5] as well as Monte 
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FIG. 2: (color-online) Correlation functions C(x,t) and C(y,t) after the quenches Qi (graphs (a) and (c)) and Q2 (graphs (b) 
and (d)) along both x and y directions. The characteristic lengths are extracted from the cutoff (see text) Co = 0.65(0.61) 
in the x(y) direction. The insets present scaling plots showing collapse of all the data for various times onto single graphs for 
each of the four case - evidence for dynamical scaling in our system. 



Carlo based renormalization group analysis [5] ■ The key 
idea in these theories is a mapping of the melting problem 
to the Kosterliz-Thoules (KT) ordering from disorder to 
quasi long-ranged order (QLRO) in the anisotropic 2d 
XY modelQSl [17]. One begins [7] by writing the elastic 
Hamiltonian, 

where u = (u x , u y ) is the displacement vector and K and 
fi are elastic constants. Note that once the external po- 
tential V(y) is applied, derivatives of u y do not appear 
since uniform translations in the transverse y direction 
costs energy. A trivial rescaling x — > \J Kx ; y ^/Jiy 
revealing the mapping to the anisotropic XY model with 
u x playing the part of the phase angle 9(r) = 2iru x /a 
and K xy = K fia 2 / Air 2 the spin-wave stiffness. Melt- 
ing in this system is governed by the unbinding of vor- 
tex like defects [T7] consisting only of dislocations pairs 
with Burgers vectors parallel to the x axis - disloca- 
tions with Burgers vectors with components in y do not 
contribute [9]. 

The kinetics of the KT transition in the 2d XY model 
following a quench from the disordered state has been 
studied in detail[TH]- The existence of dynamical scaling 
implies that the equal time correlation function, C(r, t) = 



(cos(0(O, t) — 6(r,t))) has the scaling form: 

CM) = r-v(^) (1) 

Where rj is a non-universal critical exponent (= 1/4 at 
the freezing |17j) which depends on the spin- wave stiff- 
ness and £(t) is a growing characteristic length. At large 
times £ diverges by annihilation of defect pairs - a conse- 
quence of QLRO, and C ~ r-^/(0) QZ| QI]. The rate of 
annihilation of defect pairs is controlled by free diffusion 
in 2d so that £ ~ t 1 / 2 . We show below that despite the 
mapping of the equilibrium properties of our 2d freezing 
problem to the 2dXY model, the nature of the kinetics 
of the freezing transition is distinctly different. 

To discover the kinetic processes involved in tem- 
plate assisted ordering, we carry out canonical ensemble 
Langevin dynamics simulations |19j using the fast, par- 
allelized package LAMMPS |2U|. Our system consists of 
N = 10 4 particles in a rectangular box of area A, in- 
teracting with each other at distance r through a simple 
pairwise- additive WCA interaction [21], 

U(r) = 4e [(a/r) 12 ~ (cr/r) 6 ] - e for r < r c (2) 
= otherwise 

The cutoff distance r c = 2^a and e = U(r c ). The 
units for length, energy and time are set by a, e and 
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r = (e/o^m) 1 / 2 where m is the mass. In our simulation 
all quantities are expressed in reduced units and all our 
simulations are performed at a temperature T = 0.5. 

The structure factor averaged over the simulation box 
S(qfc) = E„,meA ex PHqfc-( r m - r„)] for a 2D solid in 
a periodic potential consists of two 5 — function Bragg 
peaks at = (0, ±2n/d) and four quasi Bragg peaks at 
q/j = (±27r/a, ±27r/v / 3a), the amplitude of which, aver- 
aged over the four orientations, may be taken as a scalar 
order parameter <fi for the transition (FigJT]). We typically 
simulate the system for 4 x 10 8 steps with an integration 
time-step At = 10~ 4 with Vb = and then quench the 
system with Vq > for 7 x 10 8 steps, collecting configura- 
tions at an interval of 50000 steps. The order parameter 
(f> is averaged over the last 2000 configurations collected 
during simulation. Superimposing the last 2000 config- 
urations we observe that, for a solid phase, (f> > 0.35. 
To obtain the phase diagram shown in FigjlJ we there- 
fore took cf> — 0.35 as an (arbitrary) cut-off to obtain the 
phase boundary. This relatively crude procedure how- 
ever, is adequate for our purpose here and reproduces 
the main feature viz. re-entrant ordering in this system 
similar to that seen in earlier simulations using more so- 
phisticated met hods 019]. 

In order to obtain the appropriate correlation func- 
tion, we need the local phase 9(r). Accordingly, we 
first construct the quantities 9^ = arg{exp(iqj. • r m )} 
where q&, (fc = 1 . . . 4) are the positions of the four quasi 
Bragg peaks. The phase angles 9^ are then coarse- 
grained over blocks of size 2a x 2a centered at r, to 
obtain 9 k (r). Four correlation functions are defined as 
C k (r,t) = (cos(6> fe (0,i) - 9 k (r,t))) where the averaging 
is over the choices of the origin and over 20 independent 
quenches. Finally, C(r,t) = \J2k^k{^,t) is the average 
of these four independent functions. Since the periodic 
modulation is anisotropic, the kinetic process involved in 
freezing is expected to be different in the x and y direc- 
tions. We calculate, therefore, both C(x,t) and C(y,t). 

Our results for the correlation functions are shown 
in Fig(2] The characteristic length £ in the x and the 
y directions are defined using a suitable cutoff value 
C[£,i] = Co similar to Ref. 18J, our results do not de- 
pend crucially on the exact value of Cq. The dynamical 
scaling ansatz is seen to be valid for both quench proto- 
cols Qi and Q2 and in the x and y direction. The fitted 
value of rj(T) ss 0.1 is somewhat lower than 1/4 expected 
at freezing. This is, of course, consistent with the fact 
that 77 decreases with the stiffness of the solid which itself 
increases with the depth of quench. 

The growth of the characteristic length £ (t) is shown in 
Fig(3] This is set by the average distance between defect 
pairs[18] and shows that typically £ ~ i 1 / 4 which is the 
expected growth law if the annihilation of defect pairs 
is controlled by SFD. For the Q\ protocol at very late 
times, though, there is some indication of a crossover to 
i 1 / 2 growth especially in the transverse y direction. 
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FIG. 3: (color-online) Behaviour of the characteristic lengths 
£(i) vs t for the quench protocols Qi (a) and Q2 (b). Filled 
red (light gray) circles: £ in the x direction; open blue (dark- 
gray) circles: £ in the y direction. The straight lines show 
t 1 / 2 and i 1 / 4 behavior respectively. Note that while the time 
dependence of the Q\ quench shows some evidence for normal 
diffusion in the y direction and at late times, relaxation of 
the order parameter after the Q2 quench is always driven by 
single-file diffusion due to the strong confining effect of the 
laser potential. 



Why is the growth of the characteristic lengths set by 
SFD? The answer to this question is clear once we look 
at a typical snapshot containing a defect-antidefect pair 
as shown in Fig|4j The Burgers vectors for the two de- 
fects shown are opposite and they both lie on the same 
atomic layer. Curiously, each defect consists of a region 
of higher than average density coupled to another with 
lower density lying in an adjacent layer. The positions 
of the higher and lower densities are interchanged when 
the Burgers vector changes sign. Annihilation of defects 
therefore amounts to diffusion of particles within each 
layer making the density uniform. This is precisely the 
process involved in SFD within a channel. 

Incidentally Figj4] also explains why SFD is more 
prominent in a Q2 quench. Equilibrating the system at 
high Vq for a long time ensures that the number of parti- 
cles in each layer is nearly constant. Each layer however 
may contain large numbers of defect (anti-defect) pairs. 
However, after a quench, these defect pairs are able to 
annihilate by particles moving only within a single layer. 
Transfer of particles from one layer to another is not nec- 
essary for order to develop. On the other hand, when one 
quenches from Vq = 0, the density within each layer may 
not be the same and once defect pairs within a layer 
have annihilated with their counterparts, a few defects 
may remain where individual members of a pair exists 
on different layers. Annihilation of these defects neces- 
sarily requires particles to be transported from one layer 
to another - a process which is expected to be slow. 

In this Letter, we have described Langevin dynamics 
simulations of particles in 2d placed in a Id periodic po- 
tential. We have shown that ordering in this system oc- 
curs by diffusion of particles along the channels created 
by this potential in a single-file manner. The ordering 
length scale grows as i 1 / 4 in accordance with the predic- 
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FIG. 4: Snapshot from a close-up of our system showing the 
reduction of defect density due to pairwise annihilation of a 
dislocation with an anti-dislocation. Note that the annihila- 
tion process occurs by SFD where the order of particles within 
a layer is not altered. 



tion of SFD. How general arc our results? It is obvious 
that the nature of the template potential is crucial in 
determining the growth law. For a two dimensional po- 
tential, we expect growth to be determined by Fickean 
diffusion, though strong pinning effects may reduce the 
absolute growth rates. In actual experiments, boundaries 
may play an important role in the annihilation of defect 
pairs and may change the growth characteristics. It re- 
mains to be seen what effect, if any, open boundaries may 
have in this system which may be elucidated by future 
experiments on this system. The kinetics of ordering of 
mixtures of particlcs[22 in an external potential may also 
be interesting. 
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